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Abstract 

We will prove the equivalence of three methods, the so called energy methods, in order 
to establish the stability of an equilibrium point for a dynamical system. We will illustrate 
by examples that this result simplifies enormously the amount of computations especially 
when the stability cannot be decided with one of the three methods. 
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1 Introduction. 

Let M be a smooth manifold and 

x = f(x) (1.1) 

be a dynamical system on M given by the vector field f E X (M) and suppose x e E M is an 
equilibrium state for 1)1. lj) . i.e. f(x e ) = 0. The problem of nonlinear stability of equilibrium 
states is a very old one and the most know and remarkable results were obtained by Lyapunov 
0. They are based on finding what is called a Lyapunov function V £ C 1 (M, R) such that: 

(i) V(x e ) = 

(ii) V(x) > 0, for x ^ x e 

(iii) V < 0, where V is the derivative of V along the trajectories of 1)1.1)1 . 

In practice it is sometimes very difficult to find such a function. In many situations one can 
use constants of motion as Lyapunov functions, i.e. functions V : M — > R such that V = 0. 
This was extensively used in the context of Hamilton- Poisson systems where the Hamiltonian 
and the Casimirs of the Poisson structure are constants of motion. The methods for studying 
stability using constants of motion are the so called energy methods. The most general results 
using this methods for establishing stability can be found in !) and r6u Since in the present 
paper we are discussing local nonlinear stability we can replace, by considering a coordinate 
chart around the equilibrium x e , the manifold M with R ra , where n is the dimension of M. 

In 1965 Arnold gives the following criteria for determining nonlinear stability for an 
equilibrium point of (jl.l)) . 

Theorem 1.1 (The Arnold method [T]) Let Ci,...,Cfc <E C 2 (R n ,R) be constants of motion 
for the equation M.l]) and Fi £ C 2 (R n x R fc_1 ,R) be the smooth function given by: 

Fi(x,\i, . . . , A*, . . . , A fc ) : = Ci(x) - AiCi(x) \Ci{x) A fc C fc (x) 

where g means that the term g is omitted. If there exist constants A|, . . . , A*, . . . , A£ in R such 
that 
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(%) V x F i (x e ,\* 1 ,...,\*,...,\* k ) = 

(ii) V xx Fi(x e , A*, . . . , A*, . . . , ^* k )\ WxW is positive or negative definite, where 

k 

W := p| ker dCj(xe), 

3 = 1 

3 ± i 

then x e is nonlinear stable. 

Later, in 1985, Holm, Marsden, Ratiu and Weinstein [3] give another method for establishing 
stability of an Hamilton-Poisson system, the so called Energy-Casimir method. 

Theorem 1.2 (The Energy-Casimir method @]) L ^ Ci,C 2 , ...,C k € C 2 (R n ,R) be constants 
of motion for the equation If there exist tp%, . . . , (p%, . . . , ip k E C 2 (R, R) such that: 

(i) V x (d + (px(Ci) + ■■■ + ^d&i) + ■■■ + MCk)){x e ) = 

(ii) V 2 j. x {Ci + (pi(Ci) + • • • + tpi(Ci) + • • • + ifk{Ck)){ x e) is positive or negative definite, 
then x e is nonlinear stable. 

The above result has also an infinite dimensional analogue for Hamilton-Poison systems on 
Banach spaces, see jl]. 

Studying the stability of relative equilibria, in 1998, Ortega and Ratiu obtain, as a 
corollary of their results about stability of relative equilibria, the following theorem. 

Theorem 1.3 (The Ortega-Ratiu method [7]) Let Ct,...,C k G C 2 (lR n ,M) be constants of 
motion for the equation If there exist (pi, . . . , (pi, . . . , <p k £ C 2 (M, M) such that: 

(i) V x {d + ^i(Ci) + • • • + Sdi) + ■■■ + p k {C k ))(x e ) = 

(ii) V xx {Ci + (pi{C\) + • • • + (pi(Ci) + •••-(- (p k (C k ))(x e )\-~ is positive or negative definite, 
where 

k 

W:= f] ker(^(C,))(x e ), 

3 = 1 
3 # i 

then x e is nonlinear stable. 

The aim of our paper is to prove the equivalence of these three methods. This shows that 
when x e is an equilibrium point for (|1.1|) and we choose C±, . . . , C k as a set of constants of 
motion, if we conclude stability of x e with one of the methods, then the other two will also give 
stability of x e . Thus we can choose the most convenient method from the computational point 
of view. Since computations can become cumbersome in some examples it is important to 
know that if we cannot conclude stability of x e using the set C%, . . . , C k of constants of motion 
with one of the methods, then we cannot conclude stability of x e by applying the other two 
methods using the same set C\, . . . , C k of constants of motion. 
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2 Equivalence of the three methods 



In order for the paper to be self-contained we will start by proving Arnold's result on stability 
since in his original paper the proof was omitted. In order to do this we need the following 
preliminary results which will play a crucial role in all that follows. 

We will begin by establishing the notations and conventions to be used throughout this 
paper. A vector x £ W 1 will be considered as a column vector or a n x 1 matrix. Its transpose 
will be a row vector or a 1 x n matrix. 

Let / : M n -> M be a C 1 real valued function. The gradient of / at a point x 6 W 1 is defined 
as the column vector 

df 



V/(x) 



df 



(x) 



- dx n 

If / : R n — > M m is a vector valued map, then it will be represented as a column vector of its 
component functions fi, ■ ■ ■ , f m , namely 



fxix) 



Uf€C 



fm(x) 

l m ), then we introduce the notation 

Vf(x) := [V/i(x)...V/ m (x)], 



where V/(x) is a n x m matrix which has as columns the gradient vectors V/i(x), . . . , V/ m (x) 
Note that the transpose matrix V/(x) T is the Jacobian matrix of / at the point x G U. 



Let / : R n+k - 
following notations, 



be a C real valued function and (x, y) E 



We will use the 



V x f(x,y) 



V 2 xx f(x,y) 



V^/(x,y) 



df(x,y) ' 
dxi 

df(x,y) 
dx n 

d 2 f(x,y) 
dxidxj 

d 2 f{x,y) 



df(x,y) 



V y f(x,y) 



V|y/(s,y) 




d 2 f(x,y) 
dxidyj 



dyidyj 

For the proof of Theorem ll.ll we will need the following result that can be found in references 
[2] and • 

Proposition 2.1 Let P be a symmetric nxn matrix and Q a positive semidefinite symmetric 
n x n matrix. We assume that 

x T Px > 0, 
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for all x 6 R™, i/O satisfying x T Qx = 0. Then there exists a scalar a € R swe/t i/iai 

P + aQ > 0. 

Proof. We will prove by contradiction. Then for every integer k, there exists a vector 
Xk G R n with \\xk\\ = 1 such that: 

x\Px k + kxlQx k < 0. (2.1) 

The sequence (x^) is bounded and consequently it has a subsequence, that we will denote 
also by (xk), converging to a vector x G R™ with ||x|| = 1. Taking the limit in (|2.1j) we obtain 

x T Px+ lim (kxlQxk) < 0. (2.2) 

fc^oo 



Qx k > 0, 



Since 

the inequality (|2.2|) implies that (x^Qxk) converges to zero and hence x T Qx = 0. 
It follows from the hypothesis that x T Px > and this contradicts (|2.2[) . ■ 

Let x e be an equilibrium point for the dynamic (|l.lj) and let C±, . . . , C k G C 1 (R m ,R) be a 
set of constants of notion for the dynamic We define the following quadratic form, 

x T QjX:= ^ x T VC j (x e )(VC j (x e )) T x. (2.3) 

3 = 1 

We have the following characterization for the vector subspace W defined in Theorem ll.il 
Lemma 2.1 



T 

X 



QiX = x eW = P| ker alCj. 



3 = 1 

3 # • 



Proof. 



T 

X 



Q iX = x T X7C j (x e )(X7C j (x e )) T x = 

3 = 1 
3 # » 

^ ((VC i (x e )) T x)((VC i (x e )) T x) = 



j = i 



^ £ ((VC,(x e )fx) 2 = 



3 = 1 
3 # * 



^ (VCj(x e )) T x = 0, Vj = l,k,j^i 
44> X G If. 
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Proof of Theorem EU Let L i>ai G C 2 (R n x R k - 1 ,R) be the function denned by 
Li^iix, X\, ■ ■ ■ , Xi, . . ■ , Xk) ■= Fi(x, Ai, . . . , Aj, . . . , A&) 

+ y E [Ci(^)-Ci(^)] 2 , 

^ = i 
j ^ « 

where «j € R will be determined later. 
A simple computation shows us that 

Va;Lj )Q!i (aj, Ai, . . . , Aj, . . . , A/;) = V x Fj(aj, Ai, . . . , A*, . . . , A&) 

k 

+ ai E (C i (x)-C j (x e ))VC J (x) 

j = i 

and 

^ 7 xx-Li^cn {Xei X\ , . . . , Aj , . . . , Afc) — V xx-Fi (x e , Ai , . . . , Aj , . . . , Afc) 

+ a, ^xC j (x e )(VC j {x e )) T 

3 = 1 

3 * i 

— Pi &iQi j 

where Qi = V x Cj(x e )(VCj(x e )) T is the n x n symmetric matrix defined by (|2.3j) . 

j = i 

The hypothesis (i) implies that V x L,^ at (x e , A*, . . . , A*, . . . , X* k ) =0. As a consequence of 
the hypothesis (ii) and Proposition 12.11 we can find a* € M. such that Pi + a*Qi > and thus 
L-i,a*( x ) > for x ^ x e in a small neighborhood of the equilibrium point x e . 

Let us define now the function Vi ia * € C 2 (R n ,M) by the following relation, 

Vi,a* ( x ) = Pi,a* ( x > A*, . . . , A*, . . . , A£) — Li j(x * (x e , A*, . . . , A*, . . . , A£). 

It is easy to see that Vj )Q * is a Lyapunov function and consequently via Lyapunov's theorem 
the equilibrium state x e is nonlinear stable. ■ 

The proofs of Theorem II .21 and Theorem ll.3l can be found in the original papers [3] and [Jj. 
They are also based on finding a corresponding Lyapunov function. 
Now we will prove the main result of this paper. 

Theorem 2.2 Let C\, . . . , 6 C 2 (M n ,lR) be a set of constants of motion for the dynamic 
Then the following statements are equivalent: 

(a) hypotheses of Theorem \l.l\ hold: 

(b) hypotheses of Theorem ] 1. °A hold: 

(c) hypotheses of Theorem ] 1. 'i\ hold. 

Each of the above statements implies nonlinear stability. 
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Consider the 
arbitrary 



Proof, "(a) =4> (&)" Assume that the hypotheses of Theorem 11.11 hold 

following functions cpj : K — > R, ^ (i) = — X*t H — -(f — Cj(x e )) 2 , for j 7^ i and ay € 

* 2 _ 
for the moment, and A!- given in Theorem ll.il As in the proof of Theorem ll.il the conditions 

(i) and (ii) of Theorem 11.11 imply the conditions (i) and (ii) of Theorem 11.21 for a certain a*j 

given by Proposition 12.11 

"(b) =>■ (c)" This is obvious since positive or negative definiteness on the whole space implies 
positive or negative definiteness on the subspace W. 

"(c) (a)". Assume that the hypotheses of Theorem 11.31 hold . Let A* = —cp'j(Cj(x e )) for 
j 7^ i. It is obvious that condition (i) of Theorem 11.31 implies condition (i) of Theorem 11.11 
Also because some of Aj's might be zero we have the inclusion W C W. Then 

V 2 xx (Ci + v?i(Ci) + • • • + <fi(Ci) + ■■■ + <p k (C k ))(x e ) 

k , r 2 



- T VLC i (x e )y+ £ z T U{C l {x e )) 



1 = 1 
1 ± i 



d 2 Q(x e 



dxidxj 



y 



Yl fi{Q{x e ))z s y p 



l = 1 s,p=l 



dCi{x e ) dCi{x e ) 

dx s dx'P 



= z T vl x F i (x e ,\l,...,K,...,xt)y+ E Eyr(g t («e))^ °^ e)8 ^ a) , 

i = 1 s,p=l 
I ^ i 

for any z, y G M n . 

If we take z,y £ W the second summand will be zero and consequently condition (ii) of 
Theorem 11.31 implies condition (ii) of Theorem 11.11 ■ 

In all of the three methods the stability is decided when a certain matrix is positive or 
negative definite. Consequently, Arnold's method seems to be the most economical since it 
requires definiteness of a smaller matrix than the other two methods. 

Next we will discus the situation in which condition (i) of Theorem 11.11 is not satisfied. 
Or equivalently, when the vectors VCj(x e ), i € l,k are linear independent. Consequently, 
in a small neighborhood U Xe of x e they generate an integrable distribution whose leaves are 
the level sets of the map F := (C±, . . . ,Ck) ■ M n — > Eventually after shrinking U Xe 

all the points in U Xe are regular points for F. There exists a diffeomorphism (j) : U Xe — > 
(F^ 1 (F(x e )) n U X J x Vpr x \, where Vp( Xe ) is a small neighborhood of F(x e ) in R fc . Because 
(Ci, . . . ,Ck) are constants of motion for the dynamic (|l.lj) we obtain = (Y, 0), where 
Fel (F _1 (F(x e )) n U Xe ). If (y,z) are coordinates induced by <p on (F" 1 (F(x e ))n[/ :Ee ) xV F ( Xe ) 
from a set of coordinates around x e then the equations of motion corresponding to the vector 
field 4>*f are 



V = Y(y,z) 



0. 



(2.4) 



Moreover, <p(x e ) = (y e ,0) and y e is an equilibrium point for the dynamic generated by the 
vector field Y. The above system can be regarded as a bifurcation problem with z E Vp^ x ^ the 
bifurcation parameter. We have the following result. 
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Theorem 2.3 If (Ci, . . . , C&) are constants of motion for the dynamic 11 -1\) which are linear 
independent at the equilibrium point x e , then x e is stable for the dynamic hl.l\) if the equilibrium 
point y e is stable for the dynamic generated by the vector field Y and (y e , 0) is not a bifurcation 
point for \2.$ . 

This result was used in [3| for the stability problem of Ishii's equation. Given the conditions 
of the above theorem it is enough to study the stability of a dynamical system that has fewer 
variables. Nevertheless, the problem is not free of difficulties since one has to find a set of 
adapted coordinates around x e for the local fibration generated by the map F. 



3 Examples 

3.1 The free rigid body 

Theorem l2.2l asserts that if stability is obtained with one of the methods, then it can be obtained 
with the other two as well. Indeed, let us consider the Euler momentum equations: 



m% 



m 2 



1 


1 


h 


~h 


1 


1 


h 




1 


1 


h 





17121713 



mim-3 



m 3 = I7 — 17111712 



(M, 0, 0) is an equilibrium point and C\(mi, 1712,1713 
1 



m\ + m\ + m|) are two constants of motion. 



where I\ > I2 > I3 > 0. Then x e 
If ml m 2 2 m§\ 

2 I-7- + -j— + -7- 1 3 a nd L-2(mi, m 2 , m 3 J — - 

We study the stability of x e = (M, 0, 0), M ^ by using Arnold's method. Let F\ 
d - XC 2 , then VF x [x e ) = iff A = — . Also 



V^Fl ( X e , - 





1 1 

T 2 ~h 







1 1 



and W = Span 




1 




. It is easy to see that 



> 



. This shows that x e = (M, 0, 0), M 7^ is nonlinear stable. 

Next, we will prove the same stability result using the other two methods. We begin with 
the Energy-Casimir method. Let H v = C\ + (^(C^). The first variation is given by 

„. mi m-2 r 7713 e c , 

bHy = — — bm\ + — — dm 2 + — flm 3 + ip [miom\ + m2dm 2 + m 3 dm 3 ). 
h h h 
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Then SHJM, 0, 0) = is equivalent with ip' ( -M 



5 2 H V (M, 0,0) 



1 



h 

h~ h 



Also 



j 2 -j) (Sm 2 ) 2 • f I I ) 



+ {^M 2 )M 2 {5 mi ) 2 



is positive definite iff ip" (J^ M J > °- 

/ 1 \ 2 1 

We can take f(t) = \t M 2 I — — t and conclude that x e = (M, 0, 0), M ^ is nonlinear 

V 2 / h 

stable. 

For Ortega-Ratiu's method we can take the same constant of motion used for applying 

Arnold's method, i.e. F\ = C\ — — Ci- 

h 

3.2 Lorenz five component model 

We will show in this example that if the stability of an equilibrium point cannot be decided 
with one of the three methods then it cannot be decided with the other two either. This is 
what Theorem 12. 21 is predicting. It simplifies enormously the computations in the sense that 
if we do the computations using one of the methods and obtain that the stability cannot be 
decided, then it is useless to do the computations using the other two methods and the same 
set of constants of motion. 

To illustrate this, we will take the example of Lorenz five component model. The equations 

are 



±4 
X5 



-X2X3 + 6x2X5 

X1X3 — 6x1X5 

-XiX 2 
_X5 

h 6X1X5 

£ 



where 6, e € R*, x e 
-{x\ + 2x^ + xl + x\ + xf), and G 



(0, 0,M, 0,0), M ^ is an equilibrium point and Ci(xi, . . . , X5) 



— {x\ + x\) are constants of motion. 



We try to apply Arnold's method. Take F\ =C\ — AC2. Then Viq (x e ) = is impossible for 
any A € M. We have another possibility for choosing a constant of motion. Let F2 = C2 — \C\. 
Then VF 2 (x e ) = iff A = 0. Also 



V 2 xx F 2 (x e ,0) 



1 

10 
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and 



W = Span 



( 


1 


























1 



















i 





1 
























1 









V 

















1 


/ 



It is easy to see that S7 2 x F2(x e ,Q)\ WxW is not definite. 

Now we try to apply the Energy-Casimir method. Let Hh, = C\ + ip{C2). Then 5H^(x e ) = 
is impossible for any ip € C 2 (E, R). We take the other possibility, namely H 2 = C2 + <p(Ci). 
Then we have 

SH 2 = xiSxi + X2SX2 + (p'(x\5xi + 2x26x2 + Xz5x3 + X4SX4 + x^,5x^). 
Consequently 8H*{x e ) = iff iff {^M 2 ^j =0. Also 

5 2 Hl = (5 X1 ) 2 + (5x 2 f + v" (V) (6x 3 ? 

which is not definite. 

Finally we will try to apply Ortega- Ratiu's method. Let F = C2 + <p(Ci). We have that 
5F(x e ) = iff <f' (l;M 2 ) = and then W = E 5 . Also 



S 2 F(x e ) = 



1 
1 













ip" [ ^M 2 ) M 2 


















and consequently 5 2 F{x e )\- _ is not definite for any choice of ip G C 2 (IR, E). 
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